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Abstract. We review theoretical and experimental advances in the application 
of the perturbative QCD approach to the description of particle distributions in 
jets in the (semi) soft region with particular emphasis on HERA physics. 



1. Introduction 

In the last years hadron-jet physics has been intensively studied in e^e^, hadron- 
hadron and ep scattering processes. The advantage of HERA is that here the 
measurements can be performed at various hardness scales, from very large down 
to moderate and even small momentum transfer, thus probing the interface between 
the hard and soft physics. The data clearly show that the broad features of hadronic 
jet systems, calculated at the parton level agree surprisingly well with the measured 
ones. This proves the dominant role of the perturbative phase of jet evolution and 
supports the hypothesis of local parton-hadron duality (LPHD) jl], ||. Important 
tests of the hadroproduction dynamics come from the studies of the so-called semisoft 
region (small values of Xp = p/Ejet, but formally p ^ m/j). As long as the process of 
colour blanching and hadronization is local in the configuration space the asymptotic 
shapes of the distributions are fully predicted. The sub-asymptotic corrections are 
expressed as power series in VaJ. This is the so-called Modified Leading Logarithmic 
Approximation (MLLA) |3| which takes into account all essential ingredients of 
parton multiplication in the next-to-leading order. 

Colour coherence has led to a dramatic revision of expectations for particle 
distributions. Thus, the coherent effects in the intrajet cascades, resulting on the 
average, in the angular ordering of sequential branching, gave rise to the hump-backed 
shape of particle spectra |^ . It is not the softest particles, but those with the energies 
{Eh ^ i<;0-3-0-4^ which multiply most effectively. Due to the interjet coherence which 
is responsible for the string Q/drag Q effects in the multi-jet events, it can be verified, 
that the dynamics of the colour governs the hadroproduction. 

In this talk we focus on some selected aspects of jet physics in the (semi)soft 
region. The main goal is to illustrate some phenomcnological advances of the 
perturbative approach, for a broader outHne see, e.g.|H, ^. 
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2. Inclusive Spectra and Jet Universality 

A bright perturbative prediction is the hump-backed shape of the partonic distribution 
in the variable ^ = log ^. As a consequence of coherence soft parton multiplication 
is suppressed, and the spectrum acquires a characteristic bell-like shape ||]. It 
is amazing, that the inclusive spectra of hadrons exhibit the same shape in a large 
kinematic region in a good agreement with the MLLA-LPHD predictions. Moreover, 
the data on various hard processes demonstrate a remarkable universality assuming 
the proper choice of the cascading variables. Recall that the evolution parameter 
corresponding to the struck quark in DIS in the Breit frame is \/Q^/'2 at four- 
momentum transfer Q^. The energy scale for particle distribution in jets of energy Ej^t 
within the restricted cone Oq is Ejet&o in the small angle approximation, see The 
^-distribution can also be discussed in terms of the moments or other characteristics. 
A summary of predictions and comparison with data can be found in Here we 
discuss two recent results: the ^-distributions at HERA and the variation of the peak 
position ^* in the full energy range explored so far. 



2.1. ^ distributions at HERA 

A recent measurement of the ^ distribution by the ZEUS collaboration is shown in 
Figure |l| 0. The Gaussian shape is observed again at all jet energies Ejet = Q/2. 
The MLLA prediction is based on leading and next-to-leading order logarithmic 
approximations which become increasingly better at higher jet energies and particle 
momenta exceeding the cut-off Qo in the cascade evolution. A rather simple analytic 
expression is available in the limit when the cut-off and the QCD scale become equal, 
Qo = A ("limiting spectrum"). Because of the cut-off Qo the spectrum ends at 
^max = Y = \og{Ejf,t/Qo). The hmiting spectrum (dotted line in Figure |l|, MLLA-0) 
ends at a finite value above the peak and fits the data rather well, however, it predicts 
zero intensity beyond the limit S,max- In fact, the observed spectra in the variable 
= \og{Ejet/p) have no upper limit for p — > 0, therefore in the very soft region with 
P ^ Qo there is no agreement between theory and experiment. 

This kinematical mismatch between theory and experiment can be avoided if Qo 
is interpreted as particle mass and the theoretical prediction is taken for the parton 



energy {E = vV+Q§), or, £,e = log{Ejet/E) ^ |l3|. Then, both the theoretical 
prediction and the experimental data have the same upper boundary £,max = Y. 
Transforming back the limiting spectrum to yields the prediction MLLA-M in 
Figure |l| which fits the full spectrum quite well at the highest available energies but 
exceeds the data at lower Q. An improved but still not fully satisfactory result is 
obtained by treating the effective mass as an additional parameter (histogram in 
Figure |^). 

The discrepency of MLLA-M at large ^ can be attributed to the approximate 
form of the limiting spectrum obtained in MLLA which is not appropriate for the 
very soft particles where the coupling aaikx) becomes large. A better approximation 
for the soft particles around ^ F has been proposed already some time ago Q 
based on a modified perturbative expansion, and this analytic result has been shown 
to describe very well the ^ spectrum above the peak. These mass effects can clearly 
not be derived directly from pQCD but they are plausible and not in a contradiction 
either. Fortunately, with increasing jet energy these kinematic effects become less and 
less important as already is visible in Figure ^ 
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Figure 1. The charged particle distribution in the current fragmentation 
region in the Rieit frame for different Q bins as measured by the ZEUS 
Collaboration The full and dotted lines are the MLLA-M and MLLA-0 

predictions with A = 280 MeV respectively, the histogram is the fit of the MLLA 
with A = 280 MeV and m^f j = 230 MeV, Figure from Q . 



Similar effects appear in the determination of the moments which involve the 
integrals up to the upper boundary of the spectra. The results at the lower energies 
Q depend considerably on whether the cut-off £,max is taken into account or not. 
This leads to large discrepencies between the experimental moments |ll| and 
the predictions from the limiting spectrum 10 at low Q. On the other hand, the 
agreement is rather good if the moments are determined from the spectrum taking 
the particle mass as Qq [p^ . Again, the results from different assumptions on mass 
effects converge with increasing Q . 

Note that in spite of the known similarity between the space- and time-like 
evolution at large x, there is an essential difference between the small-x behaviour 
of the DIS structure functions and the jet fragmentation. In the former case at small 
Bjorken-a; the perturbative description at very high energies ultimately breaks down 
and the problem becomes essentially non-perturbative. There are two main limitations 
on the application of the perturbative approach: diffusion of gluons away from the hard 
scale into the infrared and unitarity constraints, see for example p5| . In the timelike 
cascades, on the contrary, due to colour coherence the density of the low Feynman-a; 
particles is not growing, and the influence of confinement dynamics reduces merely to 
an overall normalization. 
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Figure 2. Peak position of the inclusive § distribution plotted against 
di-jct mass X sinSg in comparison with the MLLA prediction (central curve); 
also shown are (in arbitrary normalization) the double logarithmic approximation 
(lower curve with asymptotic slope ~ Y/2) and expectation from cascade 
without coherence. Result by CDF Collaboration |l6| . 



Better salesmen might be tempted to claim that the observation of the hump- 
backed plateau in jet fragmentation has already clearly revealed the drastic low x- 
driven violation of DGLAP, a phenomenon which in other observables many people 
are still so desperately targeting without any success. 

2.2. Peak position S^* of distribution 

An important characteristic is the peak position The high-energy behaviour of 
this quantity is predicted [Q as 



C = Y 



C 
Y 



(1) 



with the constant term given by 



C = 



16Ncb 



= 0.2915(0.3513) 



for ' 



'■f 



= 3(5), 



with Y and Qq as given above. Furthermore a — llNc/i + 2nf/3N^, b = 
{llNc — 2nf)/3. In the large- A'^c-limit parameter C becomes independent on both 
rif and Nc and approaches its asymptotical value of C = ~ 0.23. Therefore the 
energy behaviour of ^* is determined by such a fundamental parameter of QCD as the 
celebrated ^ factor in the coefficient b. 

A nice confirmation of the perturbative picture has been presented recently by the 
CDF collaboration Q , see Figure ||. The studies of the charged particle distributions 



Multiparticle Production 



5 



were performed for a variety of di-jet masses (80 < Mjj < 630 GeV) and jet opening 
angles Qq. The shapes of the measured ^-distributions at various hardness scales 
Ejet^o turn out to be very close to the MLLA expectations. As the di-jet mass 
increases, the peak position ^* shifts towards larger values of ^ in a perfect agreement 
with the theory, whereas it joints smoothly towards the results from DIS and e^e^ 
collisions at the lower energies. 

Note that by choosing a small opening angle one can study relatively low hardness 
scales but in a more friendly environment: because of the Lorentz boost, eventually all 
particles corresponding to the "hump" become relativistic. This is a serious argument 
against the attempts to explain the particle depopulation at large ^ by the finite mass 
effects. 

3. Soft Limit of Particle Spectra 

The low momentum data provide a further confirmation of the basic ideas of QCD 
coherence and LPHD. The analysis in was performed in terms of the particle 
density for e+e^ annihilation at low momenta in quite a wide cms energy region 
(from ADONE to LEP-2). The spectra were found to be in a good agreement 
with analytical perturbative results applicable for low particle energies emphasized 
above 0, see Figure |[ The HI data [^Sj also show agreement with the 
perturbative expectations, thus confirming again the universality of soft particle 
production. 

The prediction for very low momenta p < mass is dependent to some extent on the 
treatment of kinematic mass effects as discussed in the previous section. However, the 
fact that in the soft limit p —+ the curves become flat in Figure |3[ i.e. independent 
of jet energy, is a generic result and applies for different possible implementations of 
mass effects studied in . It is a direct consequence of colour coherence: in the soft 
limit the wave length of the produced particle (gluon) is so large that it cannot resolve 
details of the primary parton jet but probes only the primary qq pair. The emission 
probability is then given by the Born term which is independent of ^/s. 

We also emphasize the observation in [0 that the approach to energy 
independence for p ^ varies with particle type. It seems that the energy 
independence is reached for momenta p < mass, namely for the heavier protons already 
at much larger p 0.6 GeV and for kaons a.t p 0.3 GeV. 

Without doubt, it is not a priori clear at all, whether one can appeal to the 
perturbative expertise when exploring the low momentum domain. However, an 
attempt to stretch the perturbative expectations to the limit of their applicability 
looks quite intriguing. This could, in principle, provide a clue for understanding 
of some conceptual problems of the LPHD. Whether or not the transition between 
two stages of jet development is smooth is a question for experiment. It is therefore 
important to test these predictions further in various ways. 

Another prediction of this kind for soft particles concerns the relative production 
in quark and gluon jets. The same arguments - dominance of Born term ~ lead one 
to expect the ratio of densities r{g/q) — -^(g — jet)/ -^{q — jet) to approach the 
limit P P 

r{9/q)^ for p-0. (2) 

This prediction would apply directly for the comparison of soft particles in and gg 
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systems. The latter system can be reahzed approximately in the process e^e~ — > qqg 
with the qq pair recoiling against the g; in this case the ratio (H) is found around 
fiQ/o) ~ 1-8 in a qualitative support of the prediction (|[), but a bit lower than 
expected, and this could be due to the non-collinearity of the three jets in the actual 
measurement. Note that the ratio of total jet multiplicities in quark and gluon jets 
stays considerably below the limit (H), both in theory and experiment, see for example 
review pO| . 

Another interesting possibility to test the prediction (|^) is in the final state of 
diffractive DIS to which we come back below. 



4. String/Drag Phenomena 

The existing data on inter-jet particle flows in e+e^ — )■ qqg events are also strongly 
in favour of soft gluon coherence. Within pQCD such effects arise from interference 
between the radiation off the q, q and g. Note that colour-related effects could play a 
valuable role as a new tool helping to distinguish the new physics signals ^. 



4.I. Particle production between jets 

Let us recall a few facts concerning comparison with the theory for the string/drag 
qqg. DELPHI pi] performed the verification for the ratio 



P _ Nqqiqqg) , . 
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of the particle densities in the interquark vaUey in the e~^e~ — s- qqg and e+e^ qq^ 
events for the F-shaped events using the g-jet tagging. The ratio in the qq angular 
interval [35°, 115°] was found to be 

i?;''P = 0.58 ± 0.06. (4) 

This is in a perfect agreement with the theoretical expectation 

K « ^2^1 « 0.61 (5) 

Another result ||2l| concerns the ratio Rg — Nqg/N^q in the q ~ g and q ~ q 
valleys in the symmetric e^e^ qqg events with the g-jet tagging. Quantitatively, 
comparing the minima at ±[50°, 70°], this ratio is found to be 

R'g'^P = 2.23 ± 0.37. (6) 

This is to be compared with the prediction [Q 

Rf « 2.4. (7) 

A clear interference pattern arises in the large- production of colourless objects, 
for instance, in V+ jet events in proton scattering (with V — j,W^ or Z). The 
hadronic antenna patterns for such processes are analogous to that in the string-drag 
case. Recently the first data on W+ jet production at the Tevatron [^ have become 
available. The interference effects are clearly seen, and they are in a good quantitative 
agreement with the perturbative calculations [Q . 

Let us emphasize that in all drag-related measurements the inter-jet flows are 
dominated by the low energy hadrons (pions with typical momenta in the few 100 
MeV range). It looks intriguing that such distant offsprings are controlled by the 
pQCD rules. 



4-. 2. Particle production perpendicular to the event plane 

An instructive test of the LPHD approach can be performed when studying the particle 
flow in the direction transverse to the 3-jet event plane, see [0. Recently DELPHI [Q 
have presented the first results on the particle yield in the transverse direction for the 
K-shaped symmetric events. Figure ^ shows the multiplicity within the cone with the 
fixed half opening angle of 30° perpendicular to the event plane. The plotted curve 
represents the lowest order formula 

Nn 1 
= 77^[2-cosei+-cosei_--^(l-cose+_)] (8) 

where the angles Qij are between q,q,g (labeled as H-,— ,1). We observe good 
agreement with the angular dependence of the perturbative prediction which provides 
a new test involving dominantly soft particles, independent of the hadronization 
models. 

The absolute normalization is also possible in terms of the radiation into a cone 
of the same size in two-jet events and this normalization is assumed in (^). It would 
correspond to three jets with 8i = in Figure |^. Changing the event topology 
from the parallel q — (qg) to anti-parallel g — (qq) configurations should increase 

the soft perpendicular radiation by the factor 9/4 as in (|^) and for 0^ = 0, 

Oi+ = Qi- = TT. In Figure ^ we are observing just the beginning of this rise as far as 
is possible in this class of symmetric events. 
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Figure 4. Multiplicity within a 30°-cone perpendicular to the event plane of 
symmetric three-jet events as a function of the angle 0i between the low energy 
jets tea]. The curve represents the perturbative prediction (H). 



It would be interesting to measure the same radiation for other three-jet topologies 
closer to the qq — g configuration such as studied by OPAL ||T^. This would test 
equation (||) closer to its maximum R± — Nc/Cp- 

Furthermore, it would be interesting to see whether the increase of radiation 
observed in Figure ^ holds down in the same way to very low particle momenta. It 
seems impossible to have a directional dependence of densities in the limit p — > 0, but 
it could become visible for momenta p ~ mass. This is the scale where the energy 
(Ejet) independence expected from the perturbative calculation establishes for tt,K 
and p particles Q . 



4-3. Conventional QCD vs. new physics mechanisms 

The antenna pattern can be used as a discriminative tool to dissect the colour structure 
of the large- events, as a way to distinguish between the conventional QCD and new 
physics mechanisms. This diagnostic power is illustrated in where the topology 
of hadronic flows corresponding to Higgs production at hadron colliders is discussed. 
There the radiation samples were compared for the signal {gg H bb -\- g) and 
background {gg bb + g) production. It was found that the main difference between 
these two came from the radiation between the final-state jets. In particular, there is, 
approximately 4/3 more radiation between these jets for the Higgs production. This 
is due to the absence of a colour connection between the quarks in the background 
process. 

Another topical example concerns central Higgs production in the events with 
double rapidity gaps, for recent discussion and references, see 
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5. Further tests with soft particles in photoproduction and DIS 

Colour coherence leads to a rich diversity of effects in jp and DIS processes (as well as 
in high-p^ events in hadronic collisions). Here we note a few lines of further studies. 



5.1. Radiation perpendicular to event plane 

The dependence of the soft perpendicular radiation on the scattering angle of the 
hard partons in various di-jet production processes have been worked out In 
particular, for small angle scattering the intensity of soft particles is larger in di-jet 
production with resolved photons (dominantly gluon exchange) by the factor Ca/Cp 
in comparison with direct photoproduction (dominant quark exchange) . Monte Carlo 
studies have shown that the expected effects should be observable under realistic 
conditions | p9| . The verification and detailed study of the angular and momentum 
dependence will provide further clues to the applicability of perturbative QCD 
calculations in this regime. 



5.2. Soft radiation in diffractive 2-jet events 

In theoretical models the diffractive cluster produced in DIS is generated through 
Pomeron-photon interactions leading to qq final state at low diffractive masses Mx 
and to qqg systems at higher masses with the gluon in the Pomeron direction and the 
qq pair in the photon direction. A recent jet analysis by the ZEUS Collaboration p^ , 
indeed, supports this view by showing the jet in the pomeron direction being broader 
than the one in photon direction. Then the soft particles in this system are expected 
again to be produced by colour octet sources and therefore the soft particle density 
is higher. This phenomenon has been observed indeed by the HI Collaboration some 



time ago |25 in the comparison of rapidity distributions of particles in the diffractive 
and the non-diffractive systems, see Figure One observes the increased central 
production of particles in the diffractive events as expected from the gluonic emitter. 
It is remarkable that the difference between the primary gluon and quark jets appears 
already at such a low mass as Mx ~ 5 GeV. The difference may become clearer if one 
restricts further to 2-jet events and to soft particles. 



5.3. Radiation in diffractive multi-jet events 

The multi-jet structure in diffractive DIS events opens up new possibilities for tests 
of coherent soft particle production. In particular, the study of diffractive three-jet 



events as in |24 allows studies complementary to the e~'"e~ case. For a given resolution 
Ucut and sufficiently high diffractive mass Mx the two-jet configuration should be 
dominated by the octet configuration like gg^ whereas the 3-jet configuration contains 
the events with a gluon jet in Pomeron direction recoiling against qq. The radiation 
into a cone perpendicular to the 3-jet plane is then governed again by the general 
formula (^ with the appropriate identification of angles. If normalized to 2-jet events 
one would replace in (H) the prefactor Nc/Cf by 1 corresponding to the replacement 
of the two-jet configuration qq by gg. 

It should be noted that particle density in the sidewise cone is not entirely energy 
{^/s = Mx) independent, this occurs only in the soft limit when the lowest order 
radiation diagrams dominate. 
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Figure 5. Rapidity distribution of charged particles in three intervals Mx of the 
difTractively produced system in DIS (HI Collaboration [p5|) in comparison with 
results from non-diffractive events at corresponding hadronic energies W = Mx ■ 



The verification of the radiation pattern (|^) in diffractive events would test both, 
the partonic interpretation of these events, and the perturbative approach to the soft 
radiation. 

6. Summary 

The existing data show that the analytical perturbative approach to multihadron 
production in QCD jets is in a remarkably healthy shape. This justifies the attempts 
to work in terms of quarks and gluons down to the low momentum scale, the region 
which is usually viewed as being entirely non-perturbative. Further tests, especially of 
the characteristic predictions from soft gluon coherence, should determine the range 
of applicability of the theory towards low momenta. 

It is clear that the agreement between parton level predictions and the data 
cannot work in all aspects. Resonances exist in nature but not in the quark gluon 
cascade. Recent studies by ZEUS Q show the qualitative agreement of perturbative 
predictions with data even for various multiparticle correlation phenomena in angular 
regions of variable size, but there is also a disagreement for correlations between 
particles at low px near the boundary Qq. In this kinematic region we encountered 
already some problems in the inclusive studies of section 2. It is therefore important 
to explore further the region of validity of this phenomenological framework. 
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The key concept of this approach is that the conversion of partons into hadrons 
occurs at low virtualities, and that it is physics of QCD branchings wliicli governs 
the gross features of the jet formation. The data demonstrate that the transition 
betwen the perturbative and non-perturbative phases is quite smooth, and that the 
spectacular coherence phenomena successfully survive the hadronization stage. 

Concluding this talk let us emphasize that, certainly, there is no mystery within 
the pQCD framework. Of real importance is that the experiment proves that the 
transformer between the perturbative and non-perturbative phases acts very smootly. 
This may (one day) shed light on the dynamics of colour confinement. 
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